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ABSTRACT 



Context. We integrate the 2D MHD ideal equations of a straight slab to simulate observational results associated with fundamental 
sausage trapped modes. 

Aims. Starting from a non-equilibrium state with a dense chromospheric layer, we analyse the evolution of the internal plasma dy- 
namics of magnetic loops, subject to line-tying boundary conditions, and with the coronal parameters described in Asai et al. (2001) 
and Melnikov et al. (2002) to investigate the onset and damping of sausage modes. 

Methods. To integrate the equations we used a high resolution shock-capturing (HRSC) method specially designed to deal appropri- 
ately with flow discontinuities. 

Results. Due to non-linearities and inhomogeneities, pure modes are difficult to sustain and always occur coupled among them so as 
to satisfy, e.g., the line-tying constraint. We found that, in one case, the resonant coupling of the sausage fundamental mode with a 
slow one results in a non-dissipative damping of the former. 

Conclusions. In scenarios of thick and dense loops, where the analytical theory predicts the existence of fundamental trapped sausage 
modes, the coupling of fast and slow quasi-periodic modes -with a node at the center of the longitudinal speed- occur contributing 
to the damping of the fast mode. If a discontinuity in the total pressure between the loop and the corona is assumed, a fundamental 
fast sausage transitory leaky regime is spontaneously produced and an external compressional Alfven wave takes away the magnetic 
energy. 

Key words. MHD; Sun: corona; Sun: oscillations; Shocks waves 



1. Introduction 

, Coronal heating remains today an open field of research. One of 
■ the most important contribution to it is attributed to active re- 
gion overdense loops which are bright in soft X-ray and EUV 
(Aschwanden et al. 2007). Over the past decade the interest in 
determining the presence of waves and oscillations that can con- 
, tribute to the heating mechanisms has increased on the basis of 
. the development of theoretical and observational studies. The 
coronal seismological remote diagnostics revealed an effective 
media to identify a wide-spectrum of fast and slow intensity 
oscillations. Measurements of characteristic periods, dispersion 
^ relations, speeds and damping times as well as different the- 
p oretical models that intend to describe non-dissipative damp- 
ing mechanisms and leakage processes of systems with high 
Reynolds number provide new remote sensing diagnostic tools 
to reveal unknown or more accurate solar physics parameters, 
e.g., magneti c field strength and transport coeffici ents (s ee, e.g. 
Aschwanden l2004l and Nakariakov and Verwichte l2005l for re- 
cent reviews). 

The pure sausage mode, modelled as a magnetic cylinder or 
a magnetic slab, is known to be a compressible MHD fast mode 
that perturb the plasma in the radial direction causing a symmet- 
ric contraction and wideni ng of t he tube without distortion of 
its axis (Edwin and Roberts 119831) . This mode was ob servation- 
ally detected in the corona by Nakariakov et al. ("2003) using the 
Nobeyama Radioheliograph. EUV and soft X-ray band periodic 
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variations of the thermal emission intensity are associated with 
density perturbations and Doppler broadening of the emission 
lines. Zaitsev and Stepanov (1989) proposed that a modulation 
of hard X-ray and white-light emission from the loop footpoints 
can occur due to the change in the loop radius and the conse- 
quent change of the mirror ratio, causing the periodic precipi- 
tation of non-thermal electrons of flaring loops. These quasi- 
periodic pulsations, triggered by sausage modes, were also as- 
sociated with flaring events of other stars (Mitra-Kraev et al. 
l2QQ5h . 

The standard theoretical models (Edwin and Roberts 119831: 

11982 : Zaitsev and Stepanov Il989h predict that the fast global 
sausage mode has two regimes separated by a cutoflT at the long 
wavenumbers: the leaky one, where the energy is lost through the 
boundaries of the loop and carried away by an external solution, 
occurring in thin and long loops when the internal Alfven speed 
is larger than the external one; and the trapped regime, occurring 
in thick and short loops, where the energy of the oscillation is 
confin ed to the interior of the loop structure. Aschwanden et al. 
(120041) showed that the sausage cutoflT imposes such a high elec- 
tron density that they can only occur in flare loops. They found 
that previously reported fast global mode observations had os- 
cillations confined to a small segment of the loop with nodes 
limiting the segment, indicating that the oscillation must be as- 
sociated with higher harmonic modes rather than fundamental 
ones. However, they found that microwave and soft X-ray ob- 
servations by Asai et al. (2001) and Melnikov et al. (2002) are 
consistent with fast sausage MHD oscillations at the fundamen- 
tal harmonic. The leaky regime has recently been studied by 
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Case A 
Case M 




Fig. 1. Divergence cleaning test. Dimensional time obtained 
multiplying by to = 1.39sec for Case A and to = 3.25sec for 
Case M. 



Pascoe et al. (l2QQ7h . who found by numerical modeling that long 
loops with sufficiently small density contrast can support global 
sausage leaky modes of detectable quality, and thus they argue 
that this mode can also be responsible of the quasi-periodic 
pulsations of flaring emission. An important issue is that the 
global sausage leaky mode period is determin ed by the Alfven 
speed outside the loop. Also, Verwichte et al. (l2009h noted that 
wave energy of loops with internal Alfven speed locally below 
the external one can tunnel through an evanescent barrier into 
the surrounding corona due to the inhomogeneity of the media 
that produces rela tive va riations of the Alfven speeds. Recently, 
Srivastava et al. (l2008l) based on the observation of multiple 
sausage oscillations suggested that fundamental post flare cool 
loops are not in equilibrium causing plasma motion along the 
loop. Pascoe et al. (l2QQ9l) investigate the efl'ects of non-uniform 
cross-sections and found that the global sausage mode is very 
efl'ectively excited in flaring loops. 

The sound (v^) and Alfven (va) speeds are associated with 
two main properties of the low corona, its compressibility and 
elasticity, respectively. The low coronal values of the parameter 
- v^/v^ « 1, imply that the sound speed is much smaller 
than the Alfven one. Due to the high Reynolds number in the 
corona and to the rapid fast mode decays, this must be firstly 
accomplished in an almost non-dissipative way. When more re- 
alistic inhomogeneous distributions of the physical quantities 
are taken into account, both in the longitudinal and radial di- 
rection (e.g., dense chromospheric layer), the pure character of 
the MHD modes is lost and the nonlinear coupling of modes 
takes place. In mediums with usual strong sudden energy release 
sources, the drain of the larger Alfvenic elastic energy to the 
lower magnetoacoustic compressive perturbations would not be 
unusual, leading to the eventual non-linear steepeni ng of m odes 
and to the formation of coronal shocks (Costa et al. l2QQ9l) . The 
quantitative evaluation of the statistical distribution and energy 
content in shock waves is important to estimate its contribution 
as a source for the coronal heating. 

In this paper our aim is to investigate the onset and damping 
of coronal sausage modes considering nonlinear and inhomoge- 
neous initial states. We simulate two studied observational flar- 
ing loops supposed capable to devel op global sausage trapped 
mo des, d escribed in Asai et al. (l2Q0lb (Case A) and Melnikov et 
al. (12002) (Case M), respectively. 



2. Numerical procedure 

The non-linear nature of MHD equations, implies that solutions 
for complex systems must be obtained by numerical means and a 
suitable numerical implementation must be constructed for this 
purpose. Even smooth initial data can easily lead to disconti- 
nuities or shock wave solutions which are not well resolved by 
standard finite difference methods and break down when the so- 
lution is not smooth enough. To simulate the plasma evolution of 
low coronal loops we used a numerical code specially designed 
to capture discontinuities via the implementation of a high reso- 
lution shock-capturing (HRSC), methods togethe r with a Runge- 
Kutta method to evolve in time (Palenzuela et al. l2009l) . 

The initial boundary conditions represent a coronal loop as a 
straight slab of enhanced density with footpoints in a dense chro- 
mospheric region and a magnetic field parallel to the axis of the 
slab. Gravity effects are neglected for simplicity and because, for 
the moment, we are not interested on curvature effects. Also the 
pressure scale height is ~ 10^ Mm. The non-dimensional ideal 
MHD equations simulated by the codes are 



dtipvi) = -dj 



dtp = djipvj) 



dte 
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e-\- p-\--B' 
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ViBiB. 
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the energy density, p is the mass density, v is the velocity, B is 
the magnetic field, p is the thermal pressure and y = 5/3 is the 
specific heat ratio. In addition there is a fictitious variable 0, cou- 
pled to the system, whose time evolution keeps under control the 
behavior of the magnetic field divergence. The line-tying of the 
magnetic field lines at the chromospheric basis (y = +L/2, L the 
length of the loop) is assured by choosing appropriate parame- 
ters that resemble the physical conditions, i.e. the change in the 
plasma parameter JS. 

The condition for the existence of a global sausage mode 
(Aschwanden et al. l2004 is 



<0.65 
oj VP, 



(1) 



OJ the width of the loop, po and pe the internal and external to the 
loop densities, respectively. The global stationary period must 
satisfy the relation: 



2L 

VAe 



2L 



2L 

vao' 



(2) 



where vao the internal Alfven speed, Vph the phase speed and, VAe 
the external Alfven speed. 

For both observational cases we used a 200^ grid of size 
xxy = 2x2, corresponding to /q = 10 Mm unities of simulation 
in Case A, and to /q = 16.3 Mm in Case M. The reference values 
used to obtain the non-dimensional equations are: = 1.39sec, 
Po = 3.56 10-^5g cm-^ = 152 G, vq = 7194km s'^ for 
Case A, and to = 3.25sec, po = 2.72 lO'^^g cm-^ = 90 G, 
vo = 5000km s"^ for Case M. The pressure reference value is 
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Time Time 

Fig. 2. Density as a function of time, across the loop width (x), and at the center of the loop length, (y = 0). Left) Case A, 
dimensional values obtained multiplying by /q = lOMm and = 1.39sec; Right) Case M, dimensional values obtained multiplying 
by /o = 16.3Mm and to = 3.25sec. 



Regions 




TaIMK] 


BaIG] 


«m[10^^ 


TmIMK] 


BmIG] 


Inner (0) 


4,3 


2.4 


152 


9.8 


9 


90 


External (e) 


0.21 


2.5 


153 


0.16 


8.8 


104 


Chromosphere (c) 


19000 


0.02 


106 


1900 


0.1 


72 



Table 1. In itial parameter values: particle number density n, te mpera ture T, magnetic field B. Left) Case A values taken from Asai 
et al. (I2QQ1I) : Right) Case M values taken from Melnikov et al. (l2QQ2h . 



calculated as po = p^v^. The loop length, the width and the inner 
density of the loops were taken in accordance with the observa- 
tional results. In the stationary stage. Case A loop has a L = 15 
Mm length and sloj = 5.6 Mm width; and Case M has a L = 24.5 
Mm length and 3. co = 5.8 Mm width. Table [T] shows the initial 
parameter values for the two trapped modes (eq.l). 

2.1. Divergence cleaning test 

As no ticed above we use a divergence cleaning scheme (Dedner 
et al. l2Q02h . We start with initial data which does not satisfy 
V • B = 0. However the departure of the calculation from this 
constraint can be evaluated through the ratio 

II V • B lb 

II B 11^1 ' 
where 

(II B Wmf = (5.,/ + B,/ + By/ + By/) dV 

Bi j indicates the derivative of Bi with respect to j. We perform 
this test for the two cases. Figure [T] shows that in five adimen- 
sional time units the evolution drives the solution to a very small 
constraint violation, < 1%. 

3. Results and discussion 

3.1. Total equilibrium pressure across the loop radius 

We first considered total equilibrium pressure between some in- 
terfaces, namely chromosphere - inner loop region and inner 



loop region - outer loop region. In the two cases, the line-tied 
loop system evolves spontaneously, in few periods (< 4) to a 
stationary regime. In the transitory, a density perturbation with 
energy draining into the exterior takes place until the sausage 
trapped regime is established. Figures|2^,b show the density per- 
turbation at the loop's center and across its transverse size (x), 
as a function of time for Case A and M, respectively. Case M 
is appreciable damped in ~ 9 periods. The sausage mode modi- 
fies the radius with a fundamental amplitude of ~ 15%7?/, (Ri the 
loop radius) in Case A and ~ ll%Ru in Case M. 

Figures [3t,b show the Fourier transform (FT) of the loop 
density for Case A and M, respectively. 
Case A 

The Case A FT (Fig.[3t) was calculated at the grid points (x, y) = 
(0.1, 0.1) and (x,y) = (0.1, 0.5), at heights of IMm from the top 
loop and 2.5Mm from the footpoint, respectively. The time is ^ = 
Tto and the dimensionless frequency is v = 1 /r. The intense peak 
of frequency 0.21 corresponds to a period of ~ 6.6sec which is 
the fundamental sausage mode described by the observations. 
The frequency 0.34 correspond to a period of 4.1 sec, associated 
to the sausage second harmonic. 

Sausage standing modes in uniform media have wavelengths 
given by the expression A = 2L/N. The departure from one of 
the ratio Pi IIP 2, between the fundamental (A/^ = 1) and second 
{N = 2) harmonic, is a measure of the longitudinal stratification 
imposed by the density contras t bet\y een the interior a nd ext erior 
of the loop (Nakariakov et al. (l2003h : McEwan et al. (l2006h ). In 
this case P1/2P2 = 0.8. 

Also, we found that there is a coupling of the sausage mode 
with slow components that are damped approximately in five pe- 
riods. This is shown in Fig. [4^, where the longitudinal velocity 
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Frequency Frequency 

Fig. 3. Density Fourier Transform for: Left) Case A, the signal is taken at the grid points (0.1,0.1) and (0.1,0.5), top loop point 
and footpoint, respectively; Right) Case M, the footpoint signal taken at the grid point (0.0, 0.5) is larger for the second peak. The 
opposite is true for the loop top signal, taken at the grid point (0.0, 0.0). 



along the loop axis is displayed. The features are a composi- 
tion of different modes, all of them with a node at the center of 
the loop length, the whole pattern resembles a second harmonic 
one indicating that this mode is predominant, as suggested by 
Nakariakov et al. (2003). The different shapes of Fig.|4^ are sim- 
ilar, so we consider the period as the time distance between two 
of these curves: ~ 17 sec. The dimensionless frequency associ- 
ated to this period is the peak 0.08, also shown in Fig. [3^. Away 
from these maxima the signal gets strongly mixed with other per- 
turbations and it is impossible to see particular features beyond 
the FT. 
CaseM 

Figure [St shows the FT spectrum of Case M for the grid points 
(x,y) = (0.0,0.0) and (x,y) = (0.0,0.5), at the top loop and at 
a height of 4.1 Mm from the footpoint, respectively. Melnikov et 
al. (2002) and Nakariakov et al. (2003) found that there are two 
main spectral components of the observed pulsations, with peri- 
ods Pi = 14 - 17sec and = 8 - Usee. Note that, as in Case A, 
the simulation of Case M shows that there are two main frequen- 
cies of characteristic periods ~ 1 5.5 sec (v = 0.21) and ~ 1 1.6sec 
(y = 0.28), respectively. Both of them are registered everywhere 
in the loop structure. However, as these authors stated, there is 
a difference in the contribution of the peaks depending on the 
location. At the legs the frequency v = 0.28 is as important as 
y = 0.21. On the other hand, at the top, the frequency v = 0.21 
has a larger amplitude than y = 0.28. Nakariakov et al. ( 20031) 
suggested that y = 0.28 corresponds to higher spatial harmonic. 
To discuss this point we performed the FT to the transverse ve- 
locity signal as an indication of transverse compression; the re- 
sult is shown in Fig.[5t. Note that, at the footpoints, the second 
peak is the most intense, the opposite occurs at the top. We found 
that the frequency 0.28 is the second sausage harmonic mode. In 
this case P1/2P2 = 0.67. Figure |4t) shows the most important 
slow mode contributions, all of them with a node at the center 
and a whole pattern resembling a second slow harmonic. There 
is a differentiated behavior in time. The first curves (dashed) 
shown in Fig. |4t) are later replaced by others (solid curves) of 
smaller amplitude and period that last not damped for at least 10 
periods. This two main periods can also be identified in Fig.[3]3 
as the frequencies y = 0.06 (54sec) and y = 0.21 (15. 5 sec). 
As in Case A, there is a slow frequency; but there is a second 
one whose frequency is coincident with the sausage fundamen- 
tal one. It seems that this slow contribution is driven by the fast 
modes of the same frequency. The damping of Case M sausage 




Fig. 4. Vy shape along the loop axis as a function of time. Top) 
Case A, dimensional quantities obtained multiplying by /q = 
lOMm, to = 1.39sec and vq = 7149Km sec"^ Bottom) Case M, 
dimensional quantities obtained multiplying by /q = 16.3Mm, 
^0 = 3.25sec and vq = 5000Km sec"^ 



signal shown in Fig.[2]3 could be due to the drain of energy from 
the sausage mode into the inner longitudinal perturbations. In 
Case A the longitudinal perturbations are not resonantly cou- 
pled to the sausage mode and thus, are not capable to damp the 
sausage mode. Figure[5]3 displays the Case M FT for the longitu- 
dinal velocity component (V3;). As a difference with Fig.[5t, due 
to the resonant slow mode contribution, the peak of the funda- 
mental sausage mode, is higher at the footpoint than at the top. 
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Fig. 5. Fourier Transform of the Case M signal of: Top) trans- 
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3.2. Plane wave mode decomposition 

To better understand the mode contribution in the observed phe- 
nomena, we perform a plane wave decomposition. This proce- 
dure gives limited information, due to the fact the wavelength of 
the perturbation is, in some cases, of the order of the background 
quantities, which make difficult to establish the real strength of 
the mode components. Some quotients between different Fourier 
components, supposedly pertaining to the same mode, have the 
expected values arising from the calculation of this standard 
mode decomposition, while others differ. However, this is de- 
pendent of the location of the signal analyzed. Thus, we com- 
pare the quotients between different field components: those re- 
sulting from the eigenvector decomposition into modes, with the 
quotients coming from the FT. The eigenvectors of the system 
were computed in previous studies (Pennisi 1993, Balsara 2001, 
Cecere et al. 2008), thus we use the ratios of the eigenvector 
components and express them as functions of the background 
quantities. As in this highly varying dynamics it is difficult to 
properly define the background quantities, we averaged the val- 
ues (in time and space) of the different fields near the center 
of the loop and near the footpoints. The comparison between 

the different mode components p/p vjp, Vy/p, BJp By/p of 3 3 equilibrium pressure across the loop radius 
Table 121 and the quotients obtamed by domg a FT of the fields: 
p/p, ^x/p, ^x/p, By/p shown in Table [3l were performed 
taking these values at the main characteristic frequencies. This 
information is extracted, in each case, from the top of the loop 
and from the footpoints. 

If for each case, we compare Table |2] with Table [3] for the 
loop top and for the footpoints, respectively, we find an overall 
reasonable agreement. However, we see that, for the loop top, 
the ratios that involve the magnetic and the velocity fields are 



Fig. 6. Density as a function of time, across the loop width (x), 
and at the center of the loop length, (y = 0). Top) Case A, 
and A;? = 0.28. Dimensional quantities obtained multiplying by 
/o = lOMm and = 1.39sec; Bottom) Case M, and Ap = 0.35. 
Dimensional quantities obtained multiplying by /q = 16.3Mm 
and ^0 = 3.25sec. 



higher than expected for y = 0.34 in Case A, and for y = 0.27 
in Case M. This is, for the second harmonic fast modes. If we 
compare the tables of the footpoints, the disagreement occurs 
only in the Melnikov slow mode (y = 0.06). On the other hand, 
the comparison between quotients of fast modes (both frequen- 
cies) shows a good accordance. These results reinforce that the 
FT analysis performed in last section, gives reliable information. 
The loop top analysis gives information about the fundamental 
sausage and the slow modes, and the footpoints analysis about 
the sausage ones. 



Post-flare loops are part of a scenario where impulsive energy 
bursts are common and become available perturbations that ex- 
cite oscillations and instabilities in magnetic structures. Thus, 
we investigated the response of the loop system to a perturba- 
tion that induces a discontinuity, i.e., an impulsive release of en- 
ergy modelled as a sudden total pressure discontinuity across 
the loop radius. We choose the dimensionless pressure pulses: 
ApA = 0.28 and ApM = 0.35. They are equivalent to a re- 
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0.02 
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0.009 
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0.0006 


0.01 
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0.01 


0.002 


0.01 
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0.27 F 


0.01 


0.005 


0.0009 
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Table 2. Main frequency strengths of the field FT quotients: X represents the FT component of the X variable. Top) Calculated near 
to the loop top and Bottom) near the footpoints. 
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Table 3. Quotients between components of the slow and fast eigenvectors: X represents the eigenvector component X, calculated 
with the averaged background values. Top) Calculated near to the loop top and Bottom) near the footpoints. 



lease of energy of ~ 500erg cm"^ and ~ 200erg cm"^ which 
are lower than the typical excess energy content in a flaring loop 
(Kirk et al, 1994). With these Rankine-Hugoniot conditions we 
obtained Figs. [6^,b. The figures show the time evolution of the 
density oscillations for Case A and Case M, respectively. Note 
that the oscillations have larger amplitudes than in the corre- 
spondent Figs.[2^,b and are rapidly damped in ~ 8 - 12 periods. 
The fundamental sausage modes have amplitudes of ~ 83%7?/, in 
Case A and ~ 100%7?/, in Case M. The damping occurs because 
the discontinuity imposes a transformation of the trapped modes 
into leaky modes. Figures [7^,b show a slice of the magnetic field 
longitudinal component where the superimposed arrows repre- 
sent the transversal velocity component for Case A and Case M, 
respectively, at the time t = 26sec and t = 44sec. As the mag- 
netic field and the velocity are frozen in quantities, the arrows are 
tracers of the energy drained away from the loop. At the center of 
the external region a compressional Alfven wave is excited and 
the perturbed magnetic energy of the loop is damped. The exci- 
tation at the footpoints is negligible due to the change in the jS 
parameter and the correspondent line-tied to the chromospheric 
bases. 



4. Conclusions 

We integrated the ideal MHD equations to simulate fundamen- 
tal trapped sausage modes observationally described (Asai et al. 



120011 (Case A) and Melnikov et al. l2002l (Case M)) incorpo- 
rating a dense chromospheric region. We could reproduce the 
observational parameters described in the literature, i.e, peri- 
ods, densities, lengths and magnetic fields. As in Melnikov et 
al. (HOOI), we found that, for both cases, there are two peak 
fast frequencies with diff'erent contributions depending on the 
location of the signal. The modes were interpreted as the fun- 
damental (more intense at the top of the loop) and the second 
sausage harmonic (more intense at the footpoints) (see Fig. [3^ 
for Case A and Fig. [Sj) and Fig. [5^ for Case M). We analyzed 
the coupling of the fundamental sausage modes with longitudi- 
nal components. In both cases we obtained slow frequency com- 
ponents, of P ~ 17sec and P ~ 54sec, respectively. Moreover, 
we found that Case M has another slow contribution, of the same 
frequency as the fundamental sausage mode, apparently driven 
by this frequency. The slow component of Case A is rapidly 
damped (~ 5 periods); whereas the -more broad slow spectra- 
of Case M lasts not damped for various periods. We suggest that 
this is due to the internal transferring of the fast mode energy into 
the slow energy which is accomplished more efficiently in the 
resonant case. Thus, we call this mechanism an internal damp- 
ing one. 

We also showed that certain initial conditions associated 
with typical flaring loop energy releases (impulsive depositions 
of energy) that could be achieved by considering jump condi- 
tions across the radio, results in the drain of energy into the exte- 
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Fig. 7. Longitudinal magnetic field component and, superim- 
posed arrows, transverse velocity component. Top) Case A, 
Bottom) Case M. 
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rior in the form of a leaky compressible Alfven mode, allowing 
the final damping of the initially trapped fundamental sausage 
mode. This external damping mechanism which is efficient to 
damp modes in an initially trapped configuration is in accor- 
dance with typical observational damping times of this modes 
lOperiods). 
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